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1.  Introduction 


Laser  velocimeters  have  found  widespread  use  for  localized  measurements  of  fluid  velocity.  The 
most  common  system  for  this  application  is  the  so-called  laser  Doppler  anemometer  (LDA).^  We 
will  consider  such  a  system,  not  for  fluid  flow  measurements,  but  for  measuring  the  velocity  of 
solid  surfaces  or  objects  large  compared  to  any  scale  defined  by  the  instrument.  In  cases  of  solid 
objects,  we  use  the  term  laser  Doppler  velocimetry  (LDV)  rather  than  anemometry.  Light  scatter¬ 
ing  in  these  cases  is  often  dominated  by  speckle  phenomena.  The  correlation  function  and  the 
corresponding  power  spectram  are  evaluated  for  different  types  of  surface  statistics  and  for  dif¬ 
ferent  parameters  defining  the  operational  mode^  of  the  instrament.  The  present  analysis  is  per¬ 
formed  using  the  method  of  generalized  ABCD  matrices.^ 

The  case  of  fully  developed  speckle  and  a  large  collection  aperture  will  give  a  power  spectram  of 
the  same  type  as  encountered  in  laser  Doppler  anemometry  with  a  large  number  of  particles.^ 

The  shape  of  the  spectram  is  then  exclusively  determined  by  the  geometry  of  the  measuring  sys¬ 
tem,  the  optical  wavelength,  and  the  velocity  of  the  scattering  object.  This  is  in  contrast  to  the  sit¬ 
uation  of  partially  developed  speckle,  or  in  cases  of  speckle  decorrelation.  In  such  cases,  the 
spectral  shape  and  the  modulation  depth  of  the  signal  will  depend  on  both  the  spatial  scales  of  the 
surface  roughness  and  on  out-of-plane  motions.  The  uncertainty  of  the  estimated  velocity  has 
been  investigated  by  Lading  and  Edwards^  in  cases  that  correspond  to  fully  developed  speckle. 
Here,  we  will  expand  the  results  so  that  partially  developed  speckle  and  speckle  decorrelation  are 
also  incorporated  in  the  analysis.  We  assume  that  photon/electron  noise  is  negligible. 

In  Sec.  2,  we  derive  general  expressions  for  the  auto  covariance  of  the  photodetector  current  in 
LDV  systems  and  present  closed-form  analytical  results  for  planar  targets  that  have  rough  sur¬ 
faces  such  that  the  reflected  optical  phase  exhibits  partial  spatial  coherence.  This  is  important  for 
measurements  on  smoother  surfaces  that  do  not  give  rise  to  fully  developed  speckle.  In  Sec.  2.2, 
we  obtain  corresponding  results  for  the  case  of  fiilly  developed  speckle  from  curved  surfaces, 
which  is  important  for  determining  the  angular  velocity  of  rotating  cylindrical  shafts.  In  Sec.  3, 
error  estimates  are  derived  that  indicate  how  well  one  can  estimate  the  velocity  of  both  planar  and 
rotating  targets  in  practice.  Finally,  in  Sec.  4,  a  comparison  with  time-of-flight  velocimetry  on 
similar  targets  is  given. 
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2.  Time-Lagged  Covariance  for  LDV  Systems 


We  consider  the  classical  LDV  system  depicted  in  Figure  1.  Two  intersecting  laser  beams  pro¬ 
duce  on  a  surface  an  interference  fringe  pattern,  which  is  imaged  onto  a  square-law  detector.  It  is 
assumed  here,  as  indicated  in  Figure  1,  that  the  x-axis  of  object-plane  coordinates  is  perpendicu¬ 
lar  to  the  interference  fringes,  and  that  the  z-axis  is  parallel  to  the  optic-axis.  It  is  also  assumed 
that  the  spacing  between  fringes  is  small  compared  to  the  overall  laser  spot  diameter  (i.e.,  the 
number  of  observable  fringes  are  much  greater  than  unity).  Furthermore,  the  diffusing  target  is 
assumed  to  move  with  a  constant  velocity,  v,  that  is  nearly  perpendicular  to  the  z-axis. 

We  seek  to  determine  the  time-lagged  (auto)  covariance  of  the  resulting  photocurrent  from  the 
detector  centered  at  the  corresponding  geometrical-image  point  of  the  center  of  the  target  spot. 
This  quantity  is  given  by 


Q(x)  =  {imt  + 1))  -  {i(t)){i(t  + 1)),  (2.1) 

where  i(t)  is  the  photocurrent  obtained  from  the  detector  at  time  t  (x  is  the  time-lag),  and  angular 
brackets  denote  the  ensemble  average  over  the  realizations  of  the  statistics  of  the  reflected  light, 
assumed  here  to  be  stationary. 


—  z 


Figure  1 .  Schematic  of  a  Laser  Doppler  Velocimeter  system.  The  transmitter 

and  receiver  are  here  shown  as  separate  units.  Often  they  are  com¬ 
bined  and  have  a  common  optical  axis. 
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The  instantaneous  photocurrent  is  given  by^ 


i(0  =  aJ*^^pW(p)/(p,0, 


(2.2) 


where  p  is  a  vector  in  the  image  plane,  /(p,t)  is  the  corresponding  image-plane  intensity  distribu¬ 
tion,  lV(p)  is  the  receiver-aperture  weighting  function,  and  a  is  a  conversion  factor  (i.e.,  power  to 
current)  given  by 


(2.3) 


where  q  is  the  electronic  charge,  T)  is  the  detector  quantum  efficiency,  v  is  the  optical  frequency, 
and  h  is  Planck’s  constant.  Here  we  model  the  (real  and  positive)  receiver  aperture  weighting 
function  by  a  circularly  symmetric  Gaussian-shaped  function  of  the  form 


lV(p)  =  exp 


(2.4) 


2 

where  Oa  is  the  1/e  radius  of  the  receiver  aperture  weighting  function. 
The  instantaneous  intensity  function  can  be  written  as 

/(P,0  = 


^/r[/^(r,/)G(r,p) 


(2.5) 


where  we  assume  that  the  detector  integration  time  is  long  compared  with  the  coherence  time  of 
the  incident  laser  light,  but  short  compared  with  the  characteristic  speckle  fluctuation  time;  Uoir,f) 
is  the  reflected  optical  field  in  the  object  plane  (assumed  here  to  be  in  the  rx-ry  plane);  and  G(r,p) 
is  the  Green’s  function  for  the  circularly  symmetric  system,  given  (to  within  an  unimportant 
phase  factor)  by 


(2.6) 
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where  A,B,  and  D  are  the  ray-matrix  components  of  the  system.^  For  the  imaging  system 
depicted  in  Figure  2,  these  quantities  are  given  by 


A  =  -f2/fi=-m,  (2.7) 

5  =  (2.8) 

ka 


D  =  -f^lf2=-l/m,  (2.9) 

2 

where  cis  the  \/e  radius  of  the  imaging  system’s  limiting  (Gaussian)  aperture,  k  is  the  optical 
wave  number,//  and/2  are  indicated  in  Figure  2,  and  m  is  the  geometrical  magnification. 
Substituting  Eqs.  (2.7)-(2.9)  into  Eq.  (2.6)  yields 


Figure  2.  Optical  diagram  for  a  LDV  system.  The  limiting  aperture  of  1/c 
radius  tris  positioned  in  the  Fourier  plane. 
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G(r,  p)  =  {-ik  /  2%B)  exp 


(x+r)2 


(2.10) 


where 


(2.11) 


(2.12) 

Equation  (2.10)  expresses  the  Green’s  function  in  terms  of  object-space  variables  w  and  x  and  is 
convenient  to  employ  in  the  ensuing  calculations. 

2.1  Planar  Targets 

We  now  assume  that  the  reflected  laser  interference  fringe  pattern,  which  is  fixed  in  space,  results 
from  reflection  off  a  partially  coherent  diffuse  planar  reflector  moving  with  uniform  velocity 
parallel  to  the  rx-axis.  In  practice,  it  is  a  good  to  assume  that  the  detector  aperture  is  larger  than 
the  imaged  spot,  and,  hence,  we  set  W(p)  =  1  in  this  section.  We  assume  that  the  amplitude  of  the 
reflection  coefficient  is  constant,  while  the  phase  exhibits  partial  spatial  coherence.  Specifically, 
we  model  the  reflected  optical  field  for  a  planar  target  as 


and 


2/i 

(0  =  -^ 
ka 


x=P. 

m 


(/^(r,0  =  (/.(r)x|/(r,r). 


(2.13) 


where  Ui(r)  is  the  incident  field.  The  mean  (diffuse)  reflected  interference  fringe  intensity  is 
given  by 


o(r)  =  Piirf  =^cos^{Kj,rx  /  2)  exp 

Tta, 


2r^ 


(2.14) 


where  Pq  is  the  reflected  power,  Og  is  the  l/e^  Gaussian  spot  radius,  and 


Kx=- 


271 


(2.15) 
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where  A,  the  fringe  period,  satisfies  the  condition  that  A  <  a®  (i.e.,  many  fringes  are  contained 
within  the  reflected  spot).  As  an  illustrative  example,  we  plot  in  Figure  3  the  reflected  intensity 
distribution  for  the  case  where  A  =  CTs/5. 

We  further  assume  that  |\|/(r,  r)|  =  p,  where  p  is  the  magnitude  of  the  reflection  coefficient.  Here 
we  assume  that  p  is  constant  and  set  it  equal  to  unity  in  the  following.  As  discussed  in  Appendix 
A,  the  correlation  function  s/\ji(r],r)\|/*(r2,/)\  is  modeled  by  a  Gaussian  function  given  by 


5y(r,,r2)  = 


f  4n^ 

^  2 

I 

1 

TT 

2 

2 

\k  y 

'•c  }) 

(2.16) 


where  is  a  measure  of  the  phase  correlation  length  of  the  target’s  surface.  Here,  we  assume 
diffuse  reflection  only.  In  particular,  we  assume  that  there  is  no  specular  component  in  the 
reflected  field  (this  is  equivalent  to  assuming  that  the  reflected  optical  phase  variance  o|  » 1  (see 
Appendix  A).  Hence,  the  reflected  field  obeys  circular  complex  Gaussian  statistics.^  In 
Appendix  B,  we  present  general  expressions  for  the  covariance  when  the  magnitude  of  the  reflec¬ 
tion  coefficient  varies  spatially.  In  the  limit  of  complete  spatial  incoherence  (i.e.,  — >  0),  we 

have  — >  (4ic/ /:^)8(r|  -r2),  where  5(r)  is  the  Dirac  delta-function.  That  is,  for  fully  developed 

speckle,  the  reflected  optical  phase  function  is  completely  random  and  delta  correlated.  For  a 


>- 

tz 

(/) 

z 

LlJ 


LU 

cr 


rx/Os 


Figure  3.  The  reflected  intensity  distribution  for  A  =  Og/S. 
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diffuser  moving  with  a  constant  velocity,  v,  we  assume  that  the  time  evolution  of  the  reflected 
phase  is  given  by 


\|/(r,  r + t)  =  \|r(r  -  vt,  t) 


(2.  17) 


That  is,  the  reflected  phase  at  any  position  in  a  coordinate  system  moving  along  with  the  diffuser 
does  not  change  in  time  (i.e.,  a  Taylor’s  hypothesis). 

On  substituting  Eqs.  (2.2)  and  (2.5)  into  Eq.  (2.1),  we  obtain  from  the  first  term  on  the  right- 
hand  side  of  Eq.  (2.1)  a  term  of  the  form 

(c/^(ri,0f/;(r2,t+T)f/*(r2,t)f/^(r;,r+T)),  (2.18) 


f  t 

where  ri,  rj,  r2,  and  r2  are  integration  variables  that  result  when  the  substitutions  alluded  to  above 
are  carried  out.  In  order  to  perform  the  statistical  average  over  the  four  indicated  reflected  fields, 
we  invoke  the  circular  complex  Gaussian  statistics  of  the  underlying  speckle  fields  to  write  the 
average  implied  by  expression  (2.18)  as 

(c/^(ri,0C/;(r2,t+'c)C/:(r2,0t/c,(r;,?+'c))  = 

(j7^(ri,r)C/;(r2,0}(c^*(ri,r + t)[/^(ri,f + 1))  (2. 1 9) 

+(f/^(ri,r)f/;(r2,t+t))(f/*(r2,0f/„(ri,r+T)}. 


From  the  discussion  preceding  Eq.  (2.17),  it  is  easily  seen  that  the  first  term  on  the  right-hand 
side  of  Eq.  (2.19)  will  yield  an  expression  that  is  identically  canceled  by  the  second  term  on  the 
right  side  of  Eq.  (2.1).  As  a  result,  we  obtain 


Q('t)  =  J  ^PiW(pi)J  dp2W{p2)K{pi,P2,'^), 


(2.20) 


where 
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(2.21) 


X  G(ri ,  Pi  )G*  (rj ,  P2  )G(ri ,  P2  )G*  (r2 ,  P2 ) 

xiy^iT^,t)UliT2,t^x)) 

x{ul(j2,t)V^{T^,t+x)). 


By  employing  Eqs.  (2.13)-(2.17),  we  obtain  that 


^(Pi.P2:'^)=  |^rif/i(ri)G(ri,pi)|jr2f/*(r2)G*(r2,P2)By(ri,r2x)  ,  (2.22) 


where 


r2x=r2+vT. 


(2.23) 


Substituting  Eqs.  (2.10),  (2.14),  (2.16),  and  (2.22)  into  Eq.  (2.20)  yields  an  integral  whose  inte¬ 
grand  contains  multiplicative  factors  of  the  form  cos(K;fri;p)cos(K;(r2tx)>  which  when  expanded, 
produce  additive  terms  of  the  form  exp[±iKj^(rijp  +  r2x  +  Vx"!)]  (  the  direct  terms),  and 
exp[±iKxVxX]  (the  cross  terms).  It  can  be  shown  that  the  direct  terms  yield,  after  oerforming  the 
integrations  over  object  space,  multiplicative  factors  of  the  form  expr-c(Kx:Cj)^r,  where  c  is  a 
dimensionless  constant  of  the  order  unity.  Because  we  assume  that  there  are  many  fringes  con¬ 
tained  within  the  reflected  spot,  =  2110^  /  A  » 1  and  hence,  the  contribution  of  the  direct 
terms  to  the  covariance  are  negligible.  On  the  other  hand,  it  can  be  shown  that  because  the  cross 
terms  are  independent  of  object  space  coordinates  they  yield  an  overall  multiplicative  factor  to 
the  covariance  that  equals  cos^  /  2)  =  [1  +  cos(Kx; v^x)]  /  2 . 

The  resulting  integration’s  have  been  performed  using  the  Mathematica  computer  program,^  with 
the  final  result  that 


p.-v  (»rfl+cos(KxVxT:) 

^  A  I  2 


(VT)^ 

“::2— j 


OJ  +  '-c 


(2.24) 


where  the  mean  current,  (i)  =  (i(t))  =  {i{t  +  x)),  is  given  by 
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{/)  =  j  W(p)J  j  ^fr2^(rhr2)i^j(riK(r2)G(ri,p)G*(r2,p) 


=  aPo 


l2/ij 


(2.25) 


and  //,  the  number  of  independent  modes  that  pass  through  the  optical  system  to  the  measure¬ 
ment  aperture,  is  given  by 


iV  = 


2Y  2^ 


V  ) 

t 

CO 


(2.26) 


Equation  (2.24)  expresses  the  LDV  covariance  of  photocurrent  as  a  function  of  the  phase  corre¬ 
lation  scale,  rc,  for  a  planar  diffuse  reflecting  target.  To  the  best  of  our  knowledge,  this  is  the  first 
analytical  model  that  includes  partial  coherence  of  the  target  for  LDV  systems. 

Examination  of  Eq.  (2.24)  reveals  that  the  effect  of  partial  coherence  of  the  target’s  surface  is  to 
increase  the  effective  size  of  the  measurement  (target)  area.  As  a  result,  it  can  be  shown  that,  for  a 
sufficiently  large  signal-to-noise  ratio,  the  number  of  observable  oscillations  in  Cj(i:)  increases 
with  increasing  r^,  hence,  increasing  the  sensitivity  of  the  system.  In  the  limit  >  0  (i.e.,  fully 
developed  speckle),  the  results  expressed  in  Eqs.  (2.24)-(2.26)  are  identical  to  that  obtained  by 
Lading  and  Edwards.^  On  the  other  hand,  for  complete  coherent  reflection  ixc  ®°),  examina¬ 
tion  of  Eq.  (2.25)  yields  that  as  expected  physically.  This  occurs  because,  for  complete 

coherent  reflection  (i.e.,  specular  reflection),  there  is  no  diffuse  component  of  the  photodetector 
current.  Thus,  for  specular  reflection,  the  LDV  covariance  of  photocurrent  is  identically  zero. 


2.2  Rotating  Targets 

In  this  section,  we  consider,  as  indicated  in  Figure  4,  a  rotating  shaft  of  radius  R  that  has  a  surface 
that  reflects  light  with  complete  incoherence  (i.e.,  rc  0).  In  all  cases  of  interest,  the  maximum 
time-lag,  T,  and  rotational  velocity  are  sufficiently  small  such  that  the  angular  separation  between 
a  point  on  the  surface  at  time  t,  and  the  corresponding  point  at  time  t  +  x  is  much  less  than  unity. 
This  is  equivalent  to  the  condition  that  the  reflected  spot  size  is  small  compared  to  the  radius  of 
the  shaft  (i.e.,  Cs/R  «1)»  ^  condition  that  is  always  met  in  practice.  For  completeness,  as  discussed 
in  the  Introduction,  we  consider  a  finite  receiver  aperture,  whose  weighting  function  is  given  by 
Eq.  (2.4).  We  than  assume  that  the  angular  velocity  of  the  shaft  is  parallel  to  the  ry-axis. 
Assuming  that  the  laser  beams  are  normally  incident  on  the  target,  the  reflected  optical  field, 

C/i(r)  in  the  initial  plane  (i.e.,  z  =  0)  at  the  coordinate  {rx,ry,0}  contains  the  multiplicative  factor 
exp(2iifcfcr/),I0  where 
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Incident 


Figure  4.  Measurement  geometry  for  a  rotating  cylindrical  shaft. 


(2.27) 


That  is,  Eq  (2.13)  becomes: 

U„(r,t)  =  f/,(r)exp(2t^fer2)vir(r,0.  (2.28) 

In  this  case,  it  can  be  shown  by  methods  similar  to  that  in  obtaining  Eq.  (2.24)  that  the  covariance 
of  the  photocurrent,  in  the  presence  of  a  rotating  shaft  and  a  finite  detector  measurement  aper¬ 
ture,  Oa  is  given  by 


(2.29) 


'  N 


'"l+cos(K;ca)^i?'t)^ 

{(OoRxf 

2  J 

exp 

where  Wq  is  the  magnitude  of  the  angular  velocity  of  the  rotating  shaft. 


_L-J_  1  (ko,(ii/Rf 

0^  (0^  +  al  tj^  +  co^ 


{i)  =  aPo 


f 

■a 

a 

1  . 

U/i. 

J 

(  .  ji 


(0+0: 


2  ^^2 


y^(a)^/m)^  +  0; 


O  J 


•-'•'a’-fe)’ 


and 


<^0 


m 


(2.30) 


(2.31) 


(2.32) 


(2.33) 


(2.34) 


(2.35) 


The  quantities  P!,  Nc,  and  So  are,  respectively,  the  number  of  independent  (optical)  modes  that 
passes  through  the  measurement  aperture  to  t^e  detector  plane,  the  co^esponding  number  of 
modes  captured  by  the  optical  system,  the  1/e^  radius  of  Ae  ^age  of  the  reflected  spot,  and  the 
detector  aperture  radius  referred  to  object-space  coordinates.  Note,  in  contrast  to  the  planar 
case,  interference  effects  due  to  the  presence  of  the  curved  target  cause  a  reduction  in  the  tempo¬ 
ral  width  of  the  covariance.  Physically,  this  arises  due  to  destractive  interference,  which  reduces 
the  effective  imaged  spot  size  and,  hence,  a  shorter  conesponding  transit  time  is  obtained. 

For  detector  apertures  that  are  large  compared  with  the  size  of  the  imaged  spot,  Eq.  (2.29) 
becomes 
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Q(t)  = 


_  (i)^  f  l  +  COs(K;(G)^/fx) 


exp 


y  J 


1  + 


(kOsOi/Ry 


1+ 


co" 


(Oa  ^  ~), 


(2.36) 


where  (i)  and  Nc  are  given  by  Eqs.  (2.31)  (with  Oa  °°)  and  (2.33),  respectively. 

As  indicated  above,  the  presence  of  the  curved  target  causes  a  reduction  in  the  temporal  width  of 
the  covariance.  For  simplicity,  consider  the  case  of  large  detector  apertures,  which  applies  for 
almost  all  systems  of  practical  concern.  Then,  examination  of  Eq.  (2.36)  reveals  that  the  effects 
of  rotation  become  significant  for  R  <  Rmim  where  Rmin  is  determined  from  the  condition  that 


ik(SX5jR^J^ 


>1. 


(2.37) 


For  example,  for  w  «  Ss,  (i.e.,  N  »  1)  Eq.  (2.37)  yields  that 


■^min  ~  kCHGr, 


while  for  fi) »  Os,  (i.e.,  IV  =  1 )  we  obtain  that 


^min  ~  kOs  ■ 


(2.38) 


(2.39) 
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3.  Error  Estimates 


Of  primary  concern  in  LDV  measurements  is  the  determination  of  the  velocity  of  a  target. 
Because  the  current  covariance  is  an  oscillating  function  of  the  time-lag  t,  velocity  information  is 
difficult  to  obtain  directly  from  the  covariance.  However,  the  corresponding  power  spectrum  is  a 
peaked  function,  from  which  velocity  information  can  be  extracted  directly  from  knowledge  of 
the  location  of  the  peak.  Therefore,  consider  the  autospectral  density  function  of  the  covariance 
of  photodetector  current,  defined  as 


oo 

S{f)= 


(3.1) 


where  Ci(x)  is  given  by  Eqs.  (2.24)  and  (2.29),  and /denotes  temporal  frequency.  Substituting 
either  Eq.  (2.24)  or  Eq.  (2.29)  into  Eq.  (3.1)  yields  a  spectram  in  the  region  of  positive  frequen¬ 
cies,  near  the  peak,  of  the  form  (excluding  the  low-frequency  region,  i.e.,  the  pedestal  term) 


U-foY 

A/^ 


(3.2) 


where 


and 


(v^  /  A  (planar  targets) 
(SioRl  A  (rotating  targets) 


A/  = 


^  v/ + 
(Of,R  /nA 


(planar  targets) 
(rotating  targets) 


(3.4) 


In  general,  in  contrast  to  planar  targets,  the  effects  of  curved  targets  results  in  an  increase  in  the 
spectral  width,  4/^  [see  Eq.  (2.30)].  Note,  in  particular,  that  if  one  knew  the  precise  location  of  the 
peak,/o,  the  corresponding  velocity  information  follows  directly  from  Eq.  (3.3).  We  next  present 
expressions  that  indicate  how  well  one  can  estimate  the  exact  location  of  the  peak  of  the  power 
spectrum  in  practice. 

For  f  0,  the  variance  of  the  estimate  of  the  autospectral  density,  Sif),  is  given  by^  ^ 


15 


(3.5) 


yar[5(/)]  =  ^(s(/)-(5(/)))^^ 

_  S\f) 

BJ  ’ 

where  T  is  the  finite  record  length,  and  Bg  is  the  resolution  bandwidth  centered  at  frequency  / 
(see  Sec.  8.5  of  Ref.  11).  The  results  given  in  Eq.  (3.5)  are  based  on  the  assumption  that  the  fil¬ 
tered  data  behave  like  bandwidth-limited  Gaussian  white  noise.  This  is  an  excellent  assumption  in 
practice  when  the  filter  resolution  bandwidth,  Be  ,  is  sufficiently  small.  The  central  limit  theorem 
applies  to  indicate  that  the  filtered  data  should  be  more  Gaussian  than  the  input  data,  and  the  fact 
that  Be  is  small  means  that  the  output  spectram  must  be  essentially  constant. 

In  order  to  obtain  an  estimate  for  the  error  in  determining  the  precise  location  of  the  peak  of 
autospectral  density  function,  we  expand  S(f)  near  f=fo,  which  yields  the  estimate 


kf)= 


if-fo9 

a/2 


Sifo)- 


Thus,  in  terms  of  estimates  of/near/o,  Var|5(/)j  is  given  by 


(3.6) 


(3.7) 


Assume  next  that  these  values  of/are  such  that/follows  a  normal  distribution  with  (/)  =  /,,  and 
variance  .  Then,  the  fourth  moment  in  Eq.  (3.7)  is  equal  to  3a^^ ,  from  which  it  follows  that 
the  variance  of  the  estimate  of  the  location  of  the  peak  in  the  autospectral  density  function  is 
given  by 


=Vi73Af2var[5(4)]^'V5(/„). 


(3.8) 


Substituting  Eq.  (3.2)  into  Eq.  (3.8)  yields  that 


(3.9) 


On  the  basis  of  the  normal  distribution  and  Eq.  (3.9),  one  can  obtain  a  confidence  interval  for  the 
unknown  true  value  of  fo,  in  terms  of  system  parameters,  based  on  any  single  estimate,  /.  Thus, 
for  example,  the  95%  confidence  interval  for  determining  the  location  of  the  peak  of  the  spec¬ 
trum  is/ =/j  ±2(5 
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4.  Discussion  and  Conclusion 


Based  on  the  results  obtained^in  this  report  for  LDV  systems  and  the  corresponding  results 
obtained  in  a  previous  article'^  treating  the  Laser-Time-of  Flight  velocimeter  (LTV),  a  comparison 
between  the  two  systems  can  be  made.  We  consider  both  results  relating  to  measurements  on  tar¬ 
gets  giving  rise  to  fully  developed  speckle  fields.  To  compare  the  two  systems,  the  transmitters 
for  the  two  systems  are  assumed  to  have  the  same  numerical  apertures,  which  make  the  number  of 
fringes  in  the  LDV  system,  Nf ,  equal  to  the  relative  distance  between  the  focused  spots  divided  by 
the  spot  radius  in  the  LTV  system.  Further,  the  numerical  apertures  of  the  detector  systems  are 
assumed  identical.  The  comparison  is  made  with  respect  to  the  error  estimates,  the  minimum 
shaft  radius,  and  the  effect  of  velocity  misalignment. 

Equation  57  of  Ref.  2  gives  the  variance  of  the  estimated  value  of  the  location  of  the  peak  of  the 
cross  covariance,  which  can  be  easily  converted  to  give  the  normalized  variance  for  estimating  the 
velocity  in  a  LTV  system: 


1 


(4.1) 


where  the  bandwidth  for  the  LTV  signal  is  determined  by  the  effective  transit  time  through  one  of 
the  illuminated  spots,  i.e.. 


Be  =  M 


1  1 


n1/2 


2  2 


(4.2) 


for  the  case  of  large  detector  apertures.  The  corresponding  normalized  variance  for  the  LDV 
system,  as  obtained  from  Eq.  (3.9),  is 


V  J  ’ 


(4.3) 


where  the  spectral  bandwidth,  for  the  case  of  large  detector  apertures,  is  also  given  by  Eq.  (4.2). 
Comparing  the  functional  dependence  of  the  normalized  variances,  Eqs.  (4.1)  and  (4.3),  reveals 
that  we  arrive  at  identical  expressions  if  the  velocity  in  the  LDV  system  is  perpendicular  to  the 
fringes.  Note  that  the  bandwidth.  Be,  is  different  for  the  two  systems.  The  LTW  system  has  a 
bandwidth  that  is  given  by  the  transit  time  through  one  of  the  focused  spots,  whereas  the  LDV 
system  has  a  corresponding  bandwidth  determined  by  the  transit  time  (which  is  a  factor  Nf  larger) 
through  the  entire  fringe  pattern.  Comparing  Eqs.  (4.1)  and  (4.3),  it  follows  that  the  variance  of 
the  estimated  velocity  is,  therefore,  a  factor  iy^l/2  smaller  for  the  LTV  system  in  comparison  to 
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the  LDV  system.  We  note  that  neither  Eq.  (4.1)  or  (4.3)  represents  fundamental  limits,5  but  even 
taking  the  detailed  signal  statistics  into  consideration  the  conclusion  appears  to  be  valid. 

Furthermore,  the  number  of  collected  modes,  Eq.  (2.33),  is  determined  by  the  spot  size  for  the 
LTV  system  and  by  the  size  of  the  fringe  pattern  for  the  LDV  system.  Identical  conditions  for 
the  detector  and  transmitter  set-ups,  therefore,  result  in  a  higher  normalized  time-lagged  auto 
covariance  for  the  LTV  system.  In  these  respects,  the  LTV  system  is  superior  to  the  LDV  system. 
A  decrease  in  resolution  is  expected  for  variations  from  optimum  alignment  for  the  LDV  system, 
whereas  the  basic  LTV  system  with  circular  spots  in  the  measuring  volume  will  lose  signal  in  this 
case. 

The  minimum  radii  of  curvature  for  the  two  systems  are  sununarized  in  Table  1.  Identical  per¬ 
formance  is  obtained  if  the  distance  between  the  spots  in  the  LTV  system  equals  the  spot  size  in 
the  LDV  system.  Both  LDV  and  LTV  systems  will  suffer  a  decrease  in  the  magnitude  of  the 
cross  covariance  and  auto  covariance  as  the  target  diameter  decreases,  while  the  temporal  widths 
of  the  respective  covariances  increase  correspondingly  with  decreasing  values  of  the  spot  size. 

Finite  spatial  coherence  of  the  reflecting  target  has  been  considered  for  the  LDV  system  only, 
showing  an  apparent  increase  [see  Eq.  (2.24)]  of  the  measuring  volume,  yielding  a  decrease  in 
the  variance  of  the  estimated  velocity.  Further,  this  is  accompanied  by  a  decrease  in  the  number 
of  optical  modes,  N,  and  therefore  an  increase  in  the  normalized  auto  covariance  and  a  slight 
decrease  in  the  magnitude  of  the  detected  signal.  Targets  giving  rise  to  fully  developed  speckle 
will,  for  both  systems,  yield  normalized  covariance  functions  that  are  inversely  proportional  to  the 
number  of  modes  collected  by  the  receiver. 

The  method  of  ABCD  matrices  has  been  applied  to  the  analysis  of  LDV  systems  in  conjunction 
with  targets  having  correlated  surface  stractures.  This  method  is  based  on  the  assumption  that  all 
apertures  are  of  Gaussian  shape.  The  resulting  Huygens-Fresnel  integrals  can  then  be  solved  ana¬ 
lytically,  whereby  later  mathematical  approximations  can  be  avoided.  The  auto  covariance  and 
power  spectrum  have  been  obtained  in  closed  analytical  form,  which  provides  a  tool  for  paramet¬ 
ric  optimization  of  the  optical  system.  System  performance  has  been  analyzed  in  that  the  vari¬ 
ance  of  the  estimated  velocity  has  been  obtained  in  terms  relating  to  the  optical  system  and  the 
target  surface  reflectance  characteristics.  Further,  measurement  restrictions  for  cylindrical  rotat¬ 
ing  targets  have  been  obtained,  and  the  results  have  been  compared  with  previous  results  for  time- 
of  flight  velocimetry. 

Table  1 .  Minimum  Radii  of  Curvature  for  Measurement  on  Cylindrical  Surfaces  in 
the  Two  Limiting  Cases  for  LDV  and  LTV 


Minimum  Radius  of  Curvature 

Time  of  Flight  Velocimetry 

Laser-Doppler  Velocimetry 

CD<Gs 

Ikdco 

kOXTs 

caxJs 

IkdGs 

kap' 
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Appendix  A 


In  this  Appendix,  we  present  a  physical  model  that  relates  the  complex  reflection  coefficient  to 
the  surface  height  fluctuations.  Following  Goodman, °  we  adopt  an  elementary  relation  that  leads 
to  tractable  results,  which  is  often  used  in  analysis  and  is  reasonably  accurate  if  the  surface  slopes 
are  small.  In  particular,  we  express  the  reflection  coefficient  as  [see  Eq.  (2.162)  and  Fig.  2.23  of 
Ref.8] 


\)r(r)  =  |\|/(r)|exp[i<t)(r)] 


(A-1) 


where. 


(l)(r)  =  A:(l  +  cosP)A(r), 


(A-2) 


h(T)  is  the  surface  fluctuation  in  surface  heights,  P  is  the  angle  between  the  direction  of  propaga¬ 
tion  of  the  laser  beams  and  the  normal  to  the  surface,  and,  for  simplicity  in  notation,  we  have  sup¬ 
pressed  the  explicit  dependence  on  time.  Because  it  is  assumed  that  /i^r)  is  a  zero-mean  wide- 
sense  stationary  random  process,  so  too  is  the  phase  angle  0(r).  If  a/j  represents  the  variance  of 
h,  then  the  variance  of  (j)  is  given  by 


o|=[/:(l+cosp)aft]^. 


(A-3) 


Assuming  that  |x|/(r)|  =  1,  and  that  the  surface  height  fluctuations  are  a  Gaussian  random  process, 
the  correlation  function  fiv|f(ri,r2)  =  ^\jr(ri)\|i*(r2)y  can  be  expressed  as 


(r,  ,r2 )  =  ^exp[j((t)(ri )  -  <l)(r2  ))]^ 

=  exp[-o^(l-iA(r,,r2))], 


(A-4) 


where  bh  is  the  normalized  correlation  function  of  the  surface  height  fluctuations. 

Equation  (A-4)  provides  us  with  a  specific  relationship  between  the  correlation  properties  of  the 
reflection  coefficient  and  the  correlation  properties  of  the  reflecting  surface.  To  proceed  further, 
we  assume  that  the  statistics  of  the  reflecting  surface  are  stationary  and  that  the  normalized  corre¬ 
lation  function  of  the  surface  heights  is  also  of  Gaussian  form. 


bh(r)  =  e\p 


(A-5) 
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where  r  =  |ri  -r2|,  and  r/j  is  the  lateral  coherence  length  of  the  surface  heights  fluctuations.  Hence, 
the  correlation  function  becomes 


5y(r)  =  exp 


-al  \-e 


(A-6) 


2 

For  almost  all  cases  of  practical  concern,  the  phase  variance  <j(|)  is  greater  than  unity  and  hence, 
can  be  expressed  as 


By^  (r)  —  exp 


(A-7) 


where 


^  =  ^  rh 

^  /:(l+cosp)a/,  ■ 


(A-8) 


is  the  phase  lateral  coherence  length. 

For  mathematical  convenience,  we  want  to  use  a  correlation  function  that  becomes  a  Dirac  delta 
function  in  the  limit  of  complete  incoherent  reflection.  That  is,  for  fully  developed  speckle,  we 
require  that  5j|,(ri,r2)  =  fl6(ri  -T2),,  where  a  is  a  constant  that  is  independent  of  the  optical  sys¬ 
tem.  The  quantity  a  can  be  determined  from  the  requirement  that  the  mean  reflected  radiance  of 
fully  developed  speckle  is  given  by  IgpAI%,  where  Iq  is  the  incident  intensity,  p  is  the  reflection 
coefficient,  and  A  is  the  area  of  the  illuminated  spot  (i.e.,  Lambert’s  law).  Hence,  it  can  be  shown 
that  a  =  47t/l^.  Based  on  these  considerations  and  Eq.  (A-7),  we  use  a  correlation  function  for 
partial  coherent  diffuse  reflection  given  by 


(  4tt'' 

2 

to 

1 

7T 

2  exp 

2  1 

) 

1) 

(A-9) 


Note,  in  the  limit  ->  0,  we  have  from  Eq.  (A-9)  that  (4jt/A:^)5(ri  -T2). 
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Appendix  B 


In  this  Appendix,  we  now  assume  that  the  magnitude  of  the  reflection  coefficient  exhibits  deter¬ 
ministic  spatial  variations  over  dimensions  of  the  reflected  laser  spot.  Specifically,  we  assume  that 

|xi/(r,0|  =  V^.  (B-1) 


where  p(r,t)  is  the  magnitude  of  the  reflection  coefficient.  Then,  following  similar  arguments  that 
led  to  Eq.  (2.22),  it  is  straightforward  to  show  that  Eq.  (2.22)  becomes 


^(Pl.P2*'t)=  Ji/rip(ri)^^^t/,(ri)G(ri,pi)|£ir2p(r2)^^^C/*(r2)G*(r2,P2)By(ri,r2x) 


(B-2) 


For  the  special  case  of  fully  developed  speckle  (i.e.,  By(ri,r2x)  =  (47i/^^j8(ri  -r2  -vx)),  Eq.  (B- 
2)  becomes 


^(Pl.P2;^)  =  ^  J  drVp(r)p(r^)  t/i(r)C/;(rOG(r,pi)G*(r„P2) 


(B-3) 


where 


=r-i-vx. 


(B-4) 
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TECHNOLOGY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for  national  security 
programs,  specializing  in  advanced  military  space  systems.  The  Corporation's  Technology 
Operations  supports  the  effective  and  timely  development  and  operation  of  national  security 
systems  through  scientific  research  and  the  application  of  advanced  technology.  Vital  to  the 
success  of  the  Corporation  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay 
aloeast  of  new  technological  developments  and  program  support  issues  associated  with  rrqpidly 
evolving  space  systems.  Contributing  c^abilities  are  provided  by  these  individu^ 
Technology  Centers: 

Electronics  Technology  Center:  Microelectronics,  VLSI  reliability,  failure 
analysis,  solid-state  device  physics,  compound  semiconductors,  radiation  effects, 
infrared  and  CCD  detector  devices,  Micro-Electro-Mechanical  Systems  (MEMS),  and 
data  storage  and  display  technologies;  lasers  and  electro-optics,  solid  state  laser  design, 
micro-optics,  optical  communications,  and  fiber  optic  sensors;  atomic  frequency 
standards,  rq)plied  laser  spectroscopy,  laser  chemistry,  atmospheric  propagation  and 
beam  control,  LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and 
evaluation,  battery  electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Center:  Evaluation  and 
characterization  of  new  materials:  metals,  alloys,  ceramics,  polymers  and  composites; 
development  and  analysis  of  advanced  materials  processing  and  deposition  techniques; 
nondestructive  evaluation,  component  failure  analysis  and  reliabUity;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and 
elevated  temperatures;  launch  vehicle  fluid  mechanics,  heat  transfer  and  flight 
dynamics;  aerothermodynamics;  chemical  and  electric  propulsion;  environmental 
chemistry;  combustion  processes;  qracecraft  structural  mechanics,  space  enviroiunent 
effects  on  materials,  hardening  and  vulnerability  assessment;  contamination,  thermal 
and  structural  control;  lubrication  and  surface  phenomena;  microengineering 
technology  and  microinstrument  development 

Space  and  Environment  Technology  Center:  Magnetospheric,  auroral  and 
cosmic  ray  physics,  wave-particle  interactions,  magnetospheric  plasma  waves; 
atmospheric  and  ionospheric  physics,  density  and  composition  of  the  upper 
atmosphere,  remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared 
astronomy,  infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  m^netosphere;  effects 
of  electrcnnagnetic  and  particulate  radiations  on  space  systems;  space  instrumentation; 
propellant  chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
atmo^heric  chemical  reactions,  atmospheric  optics,  light  scattering,  state-specific 
chemical  reactions  and  radiative  signatures  of  missile  plumes,  and  sensOT  out-of-field- 
of-view  rejectirxi. 


